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We present a status report on our calculation of the renormalization constants for the quark bilinears in 
quenched O(a) improved Wilson theory at j3 = 6.4 using quark states in Landau gauge. 



1. INTRODUCTION 

One of the leading uncertainties in current 
lattice calculations comes from the renormaliza- 
tion constants necessary to relate the lattice cur- 
rents to those in some continuum renormaliza- 
tion scheme. Until recently, the most commonly 
used method for determining these was 1-loop 
perturbation theory. In the last few years a non- 
perturbative method based on axial and vector 
Ward identities has been developed Q. This al- 
lows the determination of all scale invariant renor- 
malization constants, and all 0(a) improvement 
constants for bilinear operators. For the remain- 
ing bilinears and four fermion operators that arise 
in the weak effective Hamiltonian, the method 
of choice uses external quark states in Landau 
gauge ||. This defines the renormalization con- 
stants in the RI or MOM scheme, in which the 
value of the renormalized operators is specified 
at some fixed momentum for the external quarks. 

In this talk we present a status report of our 
calculations using quark states in Landau Gauge 
at /3 = 6.4. 

The calculation involves two quantities: the 
quark propagator, S(p) and the three point func- 
tions, C r (p, q), 
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where T is an element of the Clifford algebra, ip 
is the fermion field, and we only consider flavor 



non-singlet operators. In terms of these we define 
dSip)- 1 



= Im 



= Re 



dSjp)- 1 
dm 



= s( P )- 1 c r (p,o)s( P y 



1 R 


= V rS 


= 1 


1 R 


ry— \ry — 1pm 

= ^ 1 


= 1 


1 R 




= r 



(3) 

(4) 
(5) 

The RI scheme then corresponds to choosing a 
momentum p^ (and mass m = 0) at which the 
renormalized inverse propagator and the renor- 
malized truncated three point functions have 
their tree level values: 
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These equations, thus, define the wave function 
renormalization constant, Z^, the mass renor- 
malization constant, Z m , and the renormalization 
constants for the bilinears, Zp. The resulting 
renormalized quantities are then also the same 
those in the continuum RI scheme at renormal- 
ization scale pi 1 = p 2 . 

The most common scheme in which phe- 
nomenological results are presented is the MS, 
which can be obtained from the RI scheme re- 
sults by a perturbative calculation in the contin- 
uum. However, it is worth noting that the RI 
scheme does not respect the axial ward identities, 
and this breaking of the ward identity cannot be 
accounted for in the connection between RI and 
MS using perturbation theory Specifically, 
these identities require that 

75 (i + m fl ^-)r£ = -rs (9) 
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at q,, = m = 0, where rriR = Z m m is the renor- 
malized mass and P and A are the pseudoscalar 
and axial bilinears. Fortunately, the derivative 
terms in these equations, which contribute due to 
the Goldstone pole at zero quark mass arising out 
of chiral symmetry breaking, become irrelevant 
when the renormalization scale /x is large com- 
pared to the scale of chiral symmetry breaking ||] . 
On the other hand lattice discretization errors be- 
come large for pa larger than unity. Thus, the 
method relies on the existence of an intermediate 
window in which these two artifacts are small. 

2. DISCRETIZATION EFFECTS 

In the lattice regularization scheme operators 
receive power corrections which involve higher di- 
mension operators. This is a consequence of the 
hard cutoff, i.e. the lattice spacing a. In addition, 
because the RI scheme is defined in a fixed gauge, 
gauge-invariant and gauge-dependent operators, 
in general, mix. Application of BRS symmetry 
shows ||_that, at 0(a) and q = 0, the only correc- 
tions to ipTip are gauge invariant and of the form 

ip(WT + ryv)'0! where W is the Dirac operator 
and this correction vanishes by the equation of 
motion, ipyv = YJip = 0. Such a correction term, 
therefore, does not affect the on-shell matrix ele- 
ments of these operators, but does contribute to 
the momentum space correlators in Eq. ^| as these 
involve contributions from points where the oper- 
ator and the fermion sources overlap. 

At 0(a) only — m part of W contributes to 
the mixing; consequently, the scalar density and 
the vector current mix, as do the axial current 
and the tensor density. The pseudoscalar current 
receives only multiplicative corrections at this or- 
der. At higher orders, the violation of rotational 
and Lorentz symmetries show up, and different 
components of V^, A^, and mix amongst 
themselves. 

To take the 0(a) mixing into account, Eq. |§| is 
modified to 



where c' r is determined by the requirement that 

TrI^T'^0, (13) 

where V = S, V, T, and A for T = V, S, A, and 
T respectively. Unfortunately, as stated above, 
there is no analogous condition to determine c' P 
since 75S + S^ 1 ^ oc 75. An illustration of the 
magnitude of the mixing is given in Table [l] from 
which the c' r are determined. Similarly, at 0(a), 
the propagator gets corrections 
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where and c.^ are the gauge- variant and 

gauge-invariant corrections respectively. Meth- 
ods to determine a combination of these coeffi- 
cients have been discussed in Ref. Q and we do 
not repeat the discussion here. 

3. LATTICE PARAMETERS AND 
METHODOLOGY 

Our pilot calculation utilizes 60 lattices at f3 = 
6.4 which corresponds to a lattice scale of a -1 « 
3.85GeV || and a critical hopping parameter 
n c « 0.135796 ||. The quark propagators were 
inverted at seven values of k = 0.1280, 0.1294, 
0.1308, 0.1324, 0.1334, 0.1343, and 0.1348 using a 
clover coefficient value of csw — 1-526 JtJ. Unless 
stated otherwise, we shall only consider momenta 
such that all components are less than tt/4. 

The truncated three-point functions were cal- 
culated and projected on the Clifford basis. As 
shown in Tab. |l|, the data, in addition to the ex- 
pected terms, shows 0(a) mixing (S <-> Vi), and 
the two 0(a 2 ) effects: mixing (V3 <-> V4) and the 
difference between the terms parallel and perpen- 
dicular to the momentum. The latter effect is 
due to the violation of rotational symmetry. As 
mentioned in the previous section, we correct for 
the 0(a) mixings by determining the Cp, but, at 
pa ~ 1, the approximately 2% violation of rota- 
tional symmetry and 0.6% mixings due to higher 
order effects survive. 

To calculate the wave function renormalization, 
we Fourier transform the inverse propagator, 
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Table 1 

The various Clifford projections of the truncated 
3-point functions corresponding to the scalar 
and vector bilincars at n — 0.1294 and pa = 
(0,0,3,3) x tt/16. 



operator 



projection 



value 



v 2 
v 3 



v 4 



1 

73 
74 

7i 

72 
73 

1 

74 
74 

1 

73 



1.180(5 
0.059(3 
0.059(3 
1.067(5 
1.067(5 
1.086(3 
0.019(2 
0.006(1 
1.083(3 
0.021(2 
0.006(1 



and calculate the derivative in Eq. |^ analytically: 
T T ' = ^^2^Trj fl ^2n tl a {n} ~[[cosn u p„. (16) 

M {n} v 



We expect large 0(a 2 ) discretization errors for 
the following reason. At tree level T E = 
(1/4) C0S Pp.j which deviates from unity by 
~ 15% at pa <~ 1. To correct for this artifact, we 
divide the calculated T s by this tree-level con- 
tribution, and use the resulting T s in Eq. fj] to 
calculate Z^. Note that such a determination of 
Z^f, contains residual O(a) errors stemming from 
the use of the uncorrected propagator. 

4. RESULTS 

In Fig. |], we show the effect of improvement for 
the vector current, i.e. removing the 0(a) term by 
calculating c' v . In the absence of discretization 
errors, the calculated Zy should be independent 
of the scale \i except for non-perturbative terms. 
Non-perturbative terms are expected to fall off as 
I/// 2 at large p [Q. Instead the uncorrected data 
shows a linear rise, Zy changing by ~ 5% between 
0.8 < fia < 1.5. In contrast, the O(a) corrected 
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Figure 1. The improved (crosses) and unim- 
proved (circles) vector renormalization constants 
as a function of the renormalization scale. 



quantity]^] is constant over this range of /ia. There 
is roughly 2 — 3% uncertainty in Zy coming from 
the variation between different Lorentz compo- 
nents and different combinations of momenta in 
this range. 

Fig. U shows that the behavior of tensor renor- 
malization constant, Zt, is qualitatively similar: 
the much smaller linear rise of the unimproved 
Zt for fia > 0.8 is removed by the improvement. 
Finally, we show the behavior of the axial vec- 
tor channel in Fig. ||. We do not observe any 
improvement in Za and the data do not show a 
window in which it is independent of \xa. 

5. DISCUSSIONS 

Since we have not finished the analysis, we end 
with a few qualitative statements, (i) The O(a) 
errors induced due to mixing with the equation of 
motion operators can be corrected for by calculat- 
ing Cp. These c' T can then be compared with those 



1 arising from the use of the uncorrected propagator to 
define Zj,. 
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evaluated using the axial Ward identity |J. (ii) 
At j3 = 6.4, the presence of non-perturbative cor- 
rections forces us to work at fia ~ 0.8 or higher 
to determine the renormalization constants. At 
these large momenta, the 0(a 2 ) errors are only 
slightly smaller than the O(a) errors, and it is 
important to ascertain whether our attempts to 
correct the latter introduce unacceptably large 
0(a 2 ) errors. In particular, the correction terms 
in Eq. [l2] involve the inverse propagator, and an 
0(a) improvement in the propagator used in this 
equation changes results at 0(a 2 ). The use of 
a suitably improved propagator would make the 
0(a) correction term vanish in the chiral limit 
except for non-perturbative effects that fall off as 
1/p 2 . In practice, as already noted in Ref. ||], 
the bare inverse propagator has large 0(ap 2 ) cor- 
rections even in the chiral limit, and hence our 
improvement term does not vanish there. 

To summarize, 0(a 2 ) corrections are large and, 
at present, uncontrolled. We are investigating 
ways to improve the calculation and the results 
of this study will be reported elsewhere. 
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Figure 2. The improved (crosses) and unim- 
proved (circles) tensor renormalization constants 
as a function of the renormalization scale. 
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Figure 3. The improved (crosses) and unim- 
proved (circles) axial vector renormalization con- 
stants as a function of the renormalization scale. 



